Three Dimensional Constraint Effects on the Estimated 
ACTOD during the Numerical Simulation of Different 
Fatigue Threshold Testing Techniques 


Banavara R. Seshadri 

National Institute of Aerospace, Hampton, VA 23681 
and 

Stephen W. Smith ' 

NASA Langley Research Center, Hampton, VA 23681 

Variation in constraint through the thickness of a specimen effects the cyclic crack-tip-opening 
displacement (ACTOD). ACTOD is a valuable measure of crack growth behavior, indicating closure 
development, constraint variations and load history effects. Fatigue loading with a continual load 
reduction was used to simulate the load history associated with fatigue crack growth threshold 
measurements. The constraint effect on the estimated ACTOD is studied by carrying out three- 
dimensional elastic-plastic finite element simulations. The analysis involves numerical simulation of 
different standard fatigue threshold test schemes to determine how each test scheme affects ACTOD. The 
American Society for Testing and Materials (ASTM) prescribes standard load reduction procedures for 
threshold testing using either the constant stress ratio (R) or constant maximum stress intensity (K max ) 
methods. Different specimen types defined in the standard, namely the compact tension, C(T), and 
middle cracked tension, M(T), specimens were used in this simulation. The threshold simulations were 
conducted with different initial K max values to study its effect on estimated ACTOD. During each 
simulation, the ACTOD was estimated at every load increment during the load reduction procedure. 
Previous numerical simulation results indicate that the constant R load reduction method generates a 
plastic wake resulting in remote crack closure during unloading. Upon reloading, this remote contact 
location was observed to remain in contact well after the crack tip was fully open. The final region to 
“open” is located at the point at which the load reduction was initiated and at the free surface of the 
specimen. However, simulations carried out using the constant K max load reduction procedure did not 
indicate remote crack closure. Previous analysis results using various starting K max values and different 
load reduction rates have indicated ACTOD is independent of specimen size. A study of the effect of 
specimen thickness and geometry on the measured ACTOD for various load reduction procedures and its 
implication in the estimation of fatigue crack growth threshold values is discussed. 

Nomenclature 


«/ 

= initial crack length 

a 

= crack length 

B 

= thickness of the specimen 

C 

= stress intensity gradient, K-gradient 

da 

= crack tip element size 

E 

= Young’s modulus 

Umax 

= maximum stress intensity factor 

Rmin 

= minimum stress intensity factor 

P 

= applied pin load for C(T) specimen 

R 

= stress ratio, load ratio 

S 

= applied stress for M(T) specimen 

x,y,z 

= coordinates in x,y and z direction 

Aa 

= crack growth increment 

ACTOD 

= cyclic crack-tip-opening displacement 

AK 

= applied stress intensity factor range 
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I. Introduction 

The fatigue threshold stress intensity factor level, AK th , defines the stress intensity range, AK, at which 
a crack in metallic materials will arrest or begin to propagate. The threshold is used in the aerospace industry to 
define a durability lifetime (or safe operating time) for many components 1 . Therefore, accurate threshold data is 
critical to the safety of durability-based designs. The development of fatigue crack growth threshold data is 
standardized within organizations such as The American Society for Testing and Materials (ASTM) and 
International Standards Organization (ISO). The standards outline the experimental procedure, specimen 
geometry, and crack configurations along with tolerances on dimensions and operating parameters. The fatigue 
crack growth threshold, AK t h, for large cracks is determined by gradually reducing the stress intensity factor 
range applied during a fatigue crack growth test until a fatigue crack growth rate of 10" 10 m/cycle is achieved. 
This reducing stress intensity factor method results in a decreasing applied load as the test is performed, 
therefore, the process is known as a load reduction method. It should be noted that the rate of 10" 1() m/cycle is 
arbitrarily chosen as a value where the fatigue crack growth rate curve is nearly vertical for most materials, 
indicating a “threshold”. 

The ASTM Standard Test Method for Measurement of Fatigue Crack Growth Rates (E647) ' 
recommends that the stress, or load, ratio (R) be held constant during the required load reduction procedure to 
generate threshold data. The constant R load reduction procedure is shown schematically in Figure l.a, where 
both the maximum stress intensity factor, K max and the minimum stress intensity factor, K m j n are shown to 
decrease monotonically such that their ratio (R=K m ; n /K max ) is held constant. Experimental evidence suggests 
that the constant R load reduction procedure can result in the development of closure levels greater than the 
steady-state value 2 " 3 , which can influence the growth rate. A constant R load reduction results in larger plastic 
strains along the crack wake, created from high loads early in the test procedure, with subsequently reducing 
loads near threshold. The larger plastic strain away from the crack tip may result in contact of the plastic wake 
surfaces away from the crack tip during the unloading portion of the fatigue cycle, resulting in remote crack 
closure (Fig.l.b). In an attempt to avoid plasticity-induced crack closure along the crack wake as the threshold 
is approached, an alternative load reduction procedure has been adopted, in which the maximum stress intensity 
factor, K max , is held constant (Fig. l.a). In this technique, K max is held constant while increasing K min . thereby 
continually increasing R such that threshold measurements are made in the absence of crack closure. The 
resulting effective threshold stress intensity (AK eff ) t h is often referred to as an intrinsic measurement of fatigue 
crack growth resistance 3 " 6 . 

Two-dimensional plane stress analyses of the effects of load reduction on fatigue crack growth were 
conducted by McClung 7 , using both the finite element method 7 " 8 and a modified strip-yield model 8 . The load 
reduction procedure examined was limited to the constant R method. McClung observed elevated crack 
opening stresses as determined by the finite element analyses during the load reduction procedure. Two- 
dimensional analyses conducted by Newman 9 using a modified strip-yield model that averaged three- 
dimensional constraint effects through the thickness using an empirical constraint factor indicated that remote 
crack closure away from the crack tip can occur for a constant R load reduction procedure. This remote crack 
closure results in the crack tip opening prior to regions remote to the crack tip during loading. These elevated 
opening stress values would in turn lead to increased or non-conservative AK t h values. Daniewicz, et al. 10 
simulated both constant R and constant K max load reduction procedures using three-dimensional elastic -plastic 
finite element analyses. The three-dimensional analyses allow for a more realistic perspective of the plastic 
wake and subsequent crack opening behavior. However, the analysis was limited to large initial AK values, and 
the load reduction was confined to a few cycles due to computational limitations. Therefore, the value of AK* 
was not reached during the load reduction procedure simulation. 

The parameter, ACTOD is considered to be a fundamental parameter to characterize crack-tip 
damage 11 " 12 . To asses crack tip damage, Newman 11 evaluated ACTOD for the M(T) specimen subjected to 
various constant amplitude loading conditions. Newman considered only the plastic contribution to cyclic 
crack-tip displacement and neglected the elastic contribution leading to a variation in the results for different 


1 ASTM E647-05, Standard Test Methods for Measurement of Fatigue Crack Growth Rates 
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stress ratios. Similarly, Bichler and Pippan 12 made experimental observations of residual plastic deformation 
near the crack tip (residual CTOD) to explain the overload effects on the fatigue crack growth rate. Using 
various starting K max values and different load reduction rates, Seshadri, et al. 13-14 have shown that ACTOD is 
independent of specimen width. The estimated ACTOD values along a crack wake can be used to provide 
insight on the remote crack closure phenomenon (Fig.l.a) predicted for the constant R load reduction technique. 
Also, for the constant K max load reduction procedure, ACTOD values indicate that no closure is present 
throughout the loading. 

The main objective of this paper is to determine the constraint effects on ACTOD by numerically 
modeling fatigue crack growth in M(T) and C(T) specimens subjected to different load reduction testing 
procedures. The effect of thickness on estimated ACTOD, crack closure level and its implications on near 
threshold processes will be investigated. 

II. Finite Element Analyses 

Numerical studies have been carried out using three-dimensional elastic-plastic finite element analyses 
to model plasticity-induced crack closure under cyclic loading 8 ' 10 ' 13-15 Seshadri, et al. 16-17 have performed 
analyses with through and part-through cracks for the estimation of crack closure levels through the specimen 
thickness under constant amplitude loading, and various load reduction procedures, and have studied the effects 
of load history using three-dimensional elastic-plastic finite element analysis. Because of the very small 
element sizes required for mesh convergence, three-dimensional finite element analysis of plasticity-induced 
closure present a very computationally intensive effort. To obtain converged solutions, a minimum requirement 
is that 5-10 elements are contained within the plastic zone at any point on the crack front under the maximum 
loading 8 ' 16-17 . Crack closure and ACTOD estimates are made by monitoring nodal displacements during both the 
loading and unloading portions of cyclic loading 13-14 . The estimates are made at both “local” and “remote” 
locations. The “local” location represents a material point very close to the crack tip which characterizes the 
local crack tip region. The remote location corresponds to an initial crack tip location where earlier experimental 
evidence suggests that crack closure occurs under a constant R load reduction procedure 2-3 . Crack growth under 
cyclic loading is then simulated by loading the model with the maximum stress level of interest, and then 
releasing the nodes along the crack front to increase the crack size by one element. 

C(T) and M(T) specimens, as shown in Figure 2, exhibited two and three planes of symmetry, 
respectively, and consequently only one fourth and one eighth of the geometries were modeled using eight-node 
brick elements. Both the specimens were of 76.2 mm width and thicknesses , B of 12.7 and 2.54 mm were 
examined. The C(T) model consisted of a total of 43,752 nodes and 27,722 elements. The M(T) specimen model 
consisted of 38,712 nodes and 24,530 elements. Both specimens were analyzed using the ZIP3D 18 finite 
element code. A typical ZIP3D finite element model of the C(T) specimen is shown in Figure 3. The length of 
the crack tip elements was on the order of 10 microns. The aluminium alloy 7075-T73 was considered 
throughout the analyses, and it has a bi-linear behavior with modulus E = 71.0 GPa and flow stress ct 0 = 400 
MPa. The von Mises yield criterion and the associated flow rule were used. Small deformation theory was 
employed. The crack front was advanced one element during each cycle such that da = 10 microns. Load 
cycles were applied to simulate load reduction in the range of AK of 20 MPaVm to 1 MPaVm under both 
constant R and constant K max load reduction procedures. 

The load reduction used was defined with the following relationship: 

AK = A K o e CAa (1) 

where, Aa is the amount of crack growth during load reduction, AK 0 is the initial stress intensity factor range at 
the start of load reduction, and C (K-gradient) is a constant. For the constant R load reduction case, a value of R 
= 0.1 was used, with initial cyclic stress intensity levels, AK 0 , of 20, 15 and 10 MPaVm. For the constant K max 
load reduction case, values of K max = 22.2, 16.7 and 11.1 MPaVm were considered, with an initial value of R = 
0.1, resulting in the same initial AK values used for the constant R case. The first few load cycles were used to 
simulate constant amplitude pre-cracking with AK = AK 0 and R = 0.1. The remaining load cycles were used to 
simulate the load reduction procedures. The schematic representation of different load reduction procedures 
simulated in the current study is shown in Figure 1. The load reduction analyses presented are for a C of -500 
m -1 . While this is a greater K-gradient than prescribed for constant R load reduction in ASTM standard E 647 (- 
80 m -1 ), this rate was necessary to maintain practical computational requirements. 
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III. Results and Discussions 


The variations in local ACTOD and crack closure levels for a C(T) specimen under the constant R=0.1 
load reduction procedure with an initial AK 0 of 15 MPaVm are shown in Figure 4. The local closure and 
ACTOD values are calculated at one node, or 10 microns, behind the current crack tip location. The estimated 
ACTOD values through the thickness decrease with decreasing AK during the constant R load reduction 
procedure. At higher applied AK, ACTOD is observed to vary through the thickness as shown in Figure 4, 
where the solid black lines indicate values computed at the outer surface and dashed black lines indicate values 
computed at the mid-plane. Flowever, as the applied AK value decreases below a value of approximately 5 
MPaVm, the variation in ACTOD through the thickness is negligible. During this load reduction, local crack 
closure levels on the outer surface decrease from the steady state value of 0.38 computed at the initial AK 0 of 15 
MPaVm. For applied AK values less than 7 MPaVm, locally, the crack remains fully open through the thickness 
of the specimen (see solid and dashed red lines in Fig.4). With further reduction in applied AK to the near 
threshold regime, the estimated ACTOD value decreases while the crack remains fully open locally. 

The comparisons of remote ACTOD and crack closure level for the constant R and constant K max load 
reduction procedures are shown in Figure 5. The analysis of the C(T) specimen was performed with an initial 
AK 0 value of 15 MPaVm. In Figure 5, variations in remote ACTOD for two different through-the-thickness 
locations are provided for constant R (solid and dashed black lines in Fig. 5) and constant K max (dash-dot and 
dashed blue lines in Fig. 5) load reduction procedures. Remote crack closure levels for the constant R load 
reduction procedure (solid and dashed red lines in Fig. 5) are also provided. For the constant K max load reduction 
procedure, the crack remains fully open. No crack closure is observed for the constant R test for applied AK 
greater than 7 MPaVm; however, closure is observed to increase with decreasing AK. The initial closure is 
observed on the surface of the specimen (z = 6.35 mm) and closure at the mid-plane of the specimen is not 
observed until the applied AK approaches 4 MPaVm. While the general trend observed in Figure 5 for ACTOD 
is nearly the same for the constant R and constant K max load reduction procedures, the ACTOD values are 
greater for the K max procedure compared to the constant R load reduction data, particularly for applied AK less 
than 7 MPaVm. 

Earlier experimental 2- ' results have shown the existence of remote crack closure on the outer surface at 
higher initial AK 0 values during the constant R load reduction procedure. To understand the extent of crack 
closure at a remote location, several nodal points were chosen on the outer surface and the closure levels were 
computed during the constant R load reduction procedure. Flere, the remote location corresponds to a nodal 
point close to the initial crack tip (a;AV=0.264) just before the commencement of the load reduction procedure. 
The variations in remote crack closure levels corresponding to these remote nodal points for an M(T) specimen 
are depicted in Figure 6. The variation in crack closure levels at different remote locations versus applied AK is 
represented in Figure 6. The nodal location corresponding to the x=0.0 mm location in the figure corresponds to 
the beginning of the load reduction procedure (aj/W=0.264). The nodal points ahead of the reference point are 
indicated by positive x-coordinates, and the locations behind are indicated by negative x-coordinates. The crack 
closure level at remote locations steadily increases as the applied AK is decreased. A closure value of nearly 
one is obtained in the fatigue threshold regime; indicating the remote region (range of 200 microns) remained 
closed throughout the load cycle. From Figure 6, remote crack closure is suggested by the variation in closure 
for applied AK from approximately 19 to 10 MPaVm. As the remote closure level increases with decreasing 
AK, the local closure level decreases (solid black line in Fig. 6). For this simulation, the crack remains fully 
open locally (as indicated by the local closure level being equal to the applied R, 0.1) in the threshold regime 
and closure occurs at a remote location. Due to remote crack closure phenomenon observed under the constant 
R load reduction procedure, the local load at the crack tip will depend on the applied initial AK 0 value and load 
reduction constant, C. 

A. Effect of specimen type 

Comparison of local ACTOD and crack closure levels for M(T) and C(T) specimens under a constant 
load reduction procedure are shown in Figure 7. The comparison is made for the constant R load reduction 
procedure with an initial AK 0 of 10 MPaVm for both specimen types. Flere, local ACTOD and crack closure 
levels are shown for two different through-the-thickness locations for both the M(T) (solid and dashed black 
lines in Fig. 7) and C(T) (dash-dot and dashed blue lines in Fig. 7) specimens. Similar trends of decreasing 
ACTOD with decreasing AK and increasing ACTOD from the outer surface to the mid-plane are observed for 
both specimen types. The variation of the estimated ACTOD through the thickness is greater for the C(T) 
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specimen than is observed for the M(T) specimen, particularly at higher values of applied AK. Here, the 
ACTOD values at Z=6.35 mm and Z=0.0 mm for the M(T) specimen are bound by the values for the C(T) 
specimen at the same two locations. For the initial portion of the constant R load reduction procedure (higher 
AK), the M(T) specimen (solid and dashed red lines in Fig. 7) exhibits a higher local crack closure level through 
the thickness when compared to the C(T) (dash-dot and dashed green lines in Fig.7) specimen. As the applied 
AK becomes less than 6 MPaVm, the crack remains fully open locally for both specimen types. 


B. Effect of specimen thickness 

Numerical simulations of different load reduction procedures for both C(T) and M(T) specimens were 
carried out to understand the effects of specimen thickness on estimated ACTOD. Specimen thickness directly 
influences the calculated fatigue threshold level for various thickness specimens tested under different load 
reduction procedures 19 . The study was aimed at understanding and explaining the effects of plasticity associated 
with specimens of different thickness. Even though simulations were carried out for both C(T) and M(T) 
specimens, the results pertaining to only C(T) specimens will be discussed. Similar behavior and trends were 
observed for the M(T) specimen. A comparison of variations in local ACTOD with applied constant R load 
reduction procedure for two thicknesses are shown in Figure 8. The thicknesses examined were B = 12.7 and 
2.54 mm, representing thick and thin specimen configuration sizes, respectively. The simulations were carried 
out under identical conditions i.e., same initial applied AK 0 stress ratio R, K-gradient C and the crack tip 
element size. In general, the thin specimen is under the state of a plane stress condition compared to the thick 
specimen. This leads to an increase in plastic zone size on the outer surface of the thin specimen. In Figure 8, 
variations in ACTOD and local crack closure level during the constant R load reduction procedure simulation 
are presented for two specimen thicknesses. Due to the large plastic zone size on the outer surface (Z/B=1.0), 
the thin specimen experiences a higher crack closure level (solid and dashed red lines in Fig. 8) and decreased 
estimated ACTOD (solid and dashed black lines in Fig. 8) values compared to the thick specimen (dash-dot and 
dashed green and blue lines respectively in Fig. 8). In the interior region, (Z/B= 0.0) there is little difference 
between estimated ACTOD and local crack closure levels for thin and thick specimens. As the applied AK is 
reduced below 7 MPaVm, the crack remains fully open locally, and the estimates of ACTOD for both specimens 
are nearly identical. 

The comparisons of remote ACTOD and crack closure level for two different thickness specimens for a 
constant load reduction procedure are shown in Figure 9. The analysis of the C(T) specimen was performed with 
an initial AK 0 value of 15 MPaVm. As before, variations in remote ACTOD and closure levels are presented for 
the two different thickness specimens. The thin specimen (solid and dashed red lines in Fig. 9) experiences 
higher closure levels than the thick (dash-dot and dashed green lines in Fig. 9) specimen. Additionally, the 
estimated ACTOD is much lower for the thin specimen (solid and dashed black lines in Fig. 9) than for the thick 
specimen (dashed-dot and dashed blue lines in Fig. 9). The difference in estimated ACTOD between the thin and 
thick specimen is dependent upon the applied AK 0 value. The maximum crack closure level reaches a value of 
0.76 on the outer surface at the remote location for the thin specimen. For the thick specimen, the estimated 
maximum closure level is approximately 0.34. The higher the value of applied AK 0 , the greater the amount of 
plasticity produced for the thin specimen, which in turn leads to higher closure level at the remote location. 
Also, there is a slight variation in the estimated remote ACTOD and crack closure level through the thickness for 
both thin and thick specimens. 

The variations in closure level at several remote nodal points close to the initial crack tip just before the 
commencement of a load reduction procedure for both thin and thick specimens are shown in Figure 10. The 
nodal coordinates and their location with respect to the initial crack tip were previously described during the 
discussion of Figure 6. The nodal location corresponding to the x=0.0 mm location in the figure corresponds to 
the beginning of the load reduction procedure (a;AV=0.263). The nodal points ahead of the reference point are 
indicated by positive x-coordinates, and the locations behind are indicated by negative x-coordinates. All of 
these nodal points are located on the outer surface. For both thin and thick specimens, the remote closure level 
steadily increases as the applied AK is decreased. Due to the large plastic zone size associated with the thin 
specimens, thin specimens experience a higher closure level when compared to thick specimens as the threshold 
is approached. 

The amount of remote closure is dependent upon the applied initial AK 0 value. If the applied initial 
AK 0 value is much higher than the current value, one might expect closure levels nearly approaching one as the 
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fatigue threshold is approached, which is similar to the case shown in Figure 6. For completeness, the variation 
in local closure level with decreasing applied AK is shown for the thin (solid black line) and thick (solid green 
line) specimens in Figure 10. As the applied AK approaches the fatigue threshold regime, locally the crack tip 
remains fully open for each specimen. Flowever, at various remote locations, the closure level steadily increases 
as the applied AK is decreased in the fatigue threshold regime. 


Concluding Remarks 

Numerical estimates of ACTOD were carried out using three-dimensional elastic-plastic finite element 
analyses for both M(T) and C(T) specimens. The constant R and constant K max load reduction procedures were 
simulated to study load history effects on the fatigue crack growth threshold. Due to computational limitations, 
the analyses were performed for C = -500 m' 1 * 3 4 5 6 7 8 * 10 for all load reduction methods, which is more than 6 times greater 
than that prescribed for the constant R load reduction method. The following conclusions were made from the 
analyses. 


• The analyses confirmed the existence of remote, three-dimensional crack closure for both M(T) and 
C(T) specimens under the constant R load reduction procedure detailed in this study. 

• The analyses also indicated that the amount of remote crack closure estimated in the near threshold 
regime is primarily dependent upon the initial value of AK 0 chosen for the constant R load reduction 
procedure In contrast, under the constant K max loading, no crack closure was observed in the threshold 
regime. 

• The local ACTOD computed for the M(T) specimen are bound by the values for the C(T) specimen. 

• The effect of thickness on the estimated ACTOD and remote closure levels in the fatigue threshold 
regime were studied and a higher closure level is observed for specimens of reduced thickness. Here, 
closure is dominated by plasticity at the outer surface of the specimen. 

• The results indicate that, ACTOD is an excellent parameter to describe differences in the crack growth 
behavior observed in the laboratory based on load reduction methods and specimen configuration. 
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Figure l.a A schematic representation of different experimental threshold test techniques. 
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Figure l.b Remote crack closure phenomenon observed after constant-R load reduction procedure. 
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Figure 2. Configurations analyzed. 
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Figure 3. A typical ZIP3D finite element model of C(T) specimen. 
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Figure 4. Variation in local ACTOD with applied constant R load reduction procedure at two through- 
thickness positions (surface (6.35 mm) and mid-plane (0.0 mm)). 
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Figure 5. Comparison of remote ACTOD and closure level for different load reduction procedures for 
C(T) specimens at two through-thickness positions (surface (6.35 mm) and mid-surface (0.0 mm)). 
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Figure 6. Variation in remote crack closure level with applied constant RIoad reduction procedure for 

AK„ = 20 MPaVm 
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Material : 7075-T73 

W = 76.2 mm 
B = 12.7 mm 
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Figure 7. Comparison of variation in local ACTOD with applied constant R load reduction procedure for 
C(T) and M(T) specimens at two through-thickness positions (surface (6.35 mm) and 

mid-surface (0.0 mm)). 



1.0 


X 

0.8 | 

CO 

"o 
to 

0.6 — 
a> 
> 
a) 


a> 

i_ 

3 

(0 

o 

o 

o 

ra 

o 


0.4 


0.2 


0.0 


AK, MPaVm 


Material: 7075-T73 


C(T), W = 76.2 mm 
R = 0.1 

K max(0 f 16 66 MPa ^ m 

AK 0 = 15 MPaVm 
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Figure 8. Variation in local ACTOD and closure levels with applied constant R load reduction procedure 

for thin and thick specimens. 


11 

American Institute of Aeronautics and Astronautics 



Z/B = 1 .0 


Z/B = 0.0 


ACTOD, B=2.54 mm 



1.0 


0.8 


x 
re 

£ 

CO 

o 
CO 

0.6 

<1) 
> 
<D 


0.4 


0.2 


0.0 


re 

i_ 

3 

(A 

o 

o 

X 

o 

re 

o 


AK, MPa V m 


Material : 7075-T73 

C(T), W = 76.2 mm 
R = 0.1 

K max(o f 16 6667 MPaVm 

AK 0 = 15 MPaVm 
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Figure 9. Variation in remote ACTOD and closure levels with applied constant R load reduction 

procedure for thin and thick specimens. 
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Material : 7075-T73 

C(T), W = 76.2 mm 
R = 0.1 
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Figure 10. Comparison of variation in remote closure level on the outer surface for thin and thick 

specimens. 
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